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Abstract 

\D 

' We analyse the cosmic-ray signatures of decaying gravitino dark matter in a model- 

. independent way based on an operator analysis. Thermal leptogenesis and univer- 

sal boundary conditions at the GUT scale restrict the gravitino mass to be below 
600 GeV. Electron and positron fluxes from gravitino decays, together with the stan- 
^ ■ dard GALPROP background, cannot explain both, the PAMELA positron fraction 

and the electron -|- positron flux recently measured by Fermi LAT. For gravitino dark 
matter, the observed fluxes require astrophysical sources. The measured antiproton 
flux allows for a sizable contribution of decaying gravitinos to the gamma-ray spec- 
trum, in particular a line at an energy below 300 GeV. Future measurements of the 
gamma-ray flux will provide important constraints on possible signatures of decaying 
gravitino dark matter at the LHC. 



1 Introduction 



An unequivocal prediction of locally supersymmetric extensions of the Standard Model 
is the gravitino, the gauge fermion of supergravity [1]. Depending on the mechanism of 
supersymmetry breaking, it can be the lightest superparticle, which makes it a natural dark 
matter candidate [2]. In connection with thermal leptogenesis [3], gravitino dark matter 
has been discussed as an alternative [4] to the standard WIMP scenario [5] . 

In a class of models with small i?-parity and lepton number breaking the gravitino is no 
longer stable, but its decays into Standard Model (SM) particles are doubly suppressed by 
the Planck mass and the small i^-parity breaking parameter. Hence, its lifetime can exceed 
the age of the Universe by many orders of magnitude, and the gravitino remains a viable 
dark matter candidate [6] . Recently, it has been shown that such models yield a consistent 
cosmology incorporating nucleosynthesis, leptogenesis and gravitino dark matter [7]. 

Small /?-parity breaking can arise from spontaneous B-L breaking [7] or from left-right 
symmetry breaking [8] . Alternatively, explicit i?-parity violating couplings of heavy Ma- 
jorana neutrinos can lead to suppressed i?-parity breaking interactions in the low energy 
effective theory via the seesaw mechanism [9]. In the simplest supergravity models with 
universal gaugino masses at the Grand Unification (GUT) scale, thermal leptogenesis im- 
plies an upper bound of 600 GeV on the gravitino mass [10]. Relaxing the boundary 
conditions at the GUT scale, gravitino masses up to 1.4 TeV are possible [11]. 

Gravitino decays may lead to characteristic signatures in high-energy cosmic rays. The 
produced flux of gamma rays [6,7,12-15] and positrons [15,16] has been found to potentially 
account for the extragalactic component of the excess in the EGRET [17] and HEAT [18] 
data, respectively. Furthermore, a neutrino flux from gravitino decays is predicted [19] as 
well as a possibly observable antideuteron flux [20]. 

Recently, a steep rise in the cosmic-ray positron fraction above 10 GeV has been dis- 
covered by the PAMELA collaboration [21] whereas the observed antiproton-to-proton 
ratio [22] is consistent with previous measurements of the antiproton flux by BESS [23], 
IMAX [24] and WiZard/CAPRICE [25]. A possible explanation of this exotic positron 
source is annihilating or decaying dark matter [26], including decaying gravitinos [27,28]. 
Equally important are the recent measurements of the total electron + positron flux by 
ATIC [29], H.E.S.S. [30] and Fermi LAT [31]. 

In this paper we analyse the cosmic-ray signatures of decaying gravitino dark matter 
in a model-independent way based on an operator analysis. Consistency with the observed 
antiproton flux yields a lower bound on the gravitino lifetime. As we shall see, this de- 
termines an upper bound on the continuous gamma-ray spectrum. Following a previous 
analysis of supergravity models and leptogensis, we only consider gravitino masses below 
600 GeV. Hence, gravitino decays cannot be the cause of the anomaly observed by the 
ATIC [29] and Fermi LAT [31] collaborations. An interpretation of the PAMELA positron 
anomaly as the result of gravitino dark matter decay, on the other hand, requires gravitino 
masses above 200 GeV. 

This paper is organised as follows. In Section 2 we present a general operator analysis 
of gravitino decays and study the implications for the different branching ratios. In 
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particular we discuss the strength of the predicted monochromatic hne in the gamma-ray 
spectrum. Section 3 deals with the electron, positron and antiproton flux from gravitino 
decays and the implications for the gravitino lifetime. The results for the gamma-ray 
spectrum are discussed in Section 4, followed by our conclusions in Section 5. 



2 Gravitino decays 

i?-parity violating gravitino decays are conveniently described in terms of effective opera- 
tors. The spinor of a massive gravitino (cf. [32]) satisfies the Dirac equation 

{t^d^ - m3/2) iJu = , (2.1) 

together with the constraints 

7^^^ = , d^^^ = . (2.2) 

The mass scales multiplying the non-renormalizable operators are inverse powers of the 
Planck mass Mpi and the supersymmetry breaking gravitino mass 7713/2- This assumes for 
the masses of Standard Model particles, the gravitino mass and the masses msoft of 
other superparticles the hierarchy mgp^ -C fny2 ^ "^soft- 

Using Eqs. fl2.ll) and fl2.2p one easily verifies that the dimension-5 and dimension-6 
operators for i?-parity violating couplings of the gravitino to Standard Model particles are 
given by 

Ces = !^h'YDu4>^x + {iig'YB,, + i^gW^,) a'^^^Aj + h.c. , (2.3) 

where typically ^1,2 = 0(1/7713/2). Note, however, that in general k and the product k^i,2 
are independent parameters. For simplicity, we have suppressed the flavour indices of k, 
/ and ^1^2- The covariant derivative involves the U(l) and SU(2) gauge fields and W^, 
respectively, 

D^ = d, + ig'YB^ + tgW,, ^[0] = , 1^^ = Vw^^^ (2.4) 
with the corresponding field strengths 

B,, = d,B, - d,B^ , w,, = d,w, - d,w, + i[w,, w,] . (2.5) 

In the unitary gauge, the Higgs and lepton doublets read 

*1 . (2.6) 



/ v^. 

where /i, p and e denote Higgs boson, neutrino and charged lepton, respectively. 
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From Eq. (12.3^ one easily obtains the couplings of leptons and gravitino to Higgs and 
gauge bosons, which are responsible for the two-body gravitino decays, 



V 2Mpi 

+ iV2mwiwd^M;;ei''(y^''i^x} + h.c. , (2.7) 



with 
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iz = sin + cos 6'w6 , = 

g' 

= sm 6'w cos 9w (6 - 6) , sm 9w = , ,^ ■ (2.8) 

V9 +9 

Note that the gauge boson couplings satisfy the relation 

^z + tanew^^ = ^w ■ (2.9) 

For ^1 = one has = and ^z = ^w- 

The interaction Lagrangian fl2.7p coincides with the one obtained from bilinear i?-parity 
breaking [19,33] if parameters are properly matchecQ. Using the results of [19,33] we then 
obtain the partial gravitino decay widths 

r(V'3/2 ^ hvi) = Pt , (2.10) 

^Pl V ^3/2 



r(^3/2 - ^i^.) = 3g^^2f2 t + 16^— + 6mi|ez.|'Fz ) , (2.12) 



F V^3/2 ^ e+ = 2 + 16 Giy + 6m;^^iyi i^w 

1927rMpj \ m3/2 



(2.13) 



where the subscript i denotes the generation index. The functions Ha, Ga and Fa 
{a = h, Z, W) are given by [19] 

/3a = 1 - ^ , (2.14) 

""^3/2 

if, = 1 + 10^ + , (2.15) 

"^3/2 "^3/2 



^ The relations are Ki = {vi)/v, U^^ ~ —mz^z, U^y^ = —mwS.w, ~ ~itiz£,^, C/^^^^i'^/^ 
U^^^cosP + mlJiml^-ml) = 1. 



1 M2 

2 777,3^2 
2 M2 1 



G'a = l + on^' (2-16) 



F„ = l + -^ + -^. (2.17) 

J 777,3^2 'J "^3/2 

As expected, one has for maC,a ~ '^a/"^3/2 ^ 1, 

r(v^3/2 ^ /iz/.) ^ r(v^3/2 ^ ^i^.) ^ ^r(v^3/2 ^ w^^e^) . (2.18) 

The decay width T {1^3/2 — > 7z/) is of particular interest since it determines the strength 
of the gamma hne at the end of the continuous spectrum. As discussed above, this decay 
width is model-dependent. Contrary to the continuous part of the spectrum it can vanish, 
which is the case for ^1 — ^2 = 0. Generically, without such a cancellation, one obtains for 
the branching ratio using Eqs. ([23]) and (IXTOD - (|2A3|) : 

BR(V'3/2 ^ 1^) ~ 0.3 sin^ ew ( ^) ' ~ 0.02 ( ^''^ ) ' . (2.19) 

V"^3/2/ V "^3/2 / 

This estimate will be used in Section 4 where the gamma-ray spectrum is discussed. 

Another phenomenologically important issue is the flavour structure of gravitino de- 
cays, i.e., the dependence of the parameters ^^j, ^zi and on the generation index. In 
models of bilinear i?-parity breaking, for instance, this information is encoded in the mixing 
parameters /ij, 

AC = i2,HuLi , (2.20) 

where Hu and are Higgs and lepton-doublet superfields, respectively. In models of flavour 
the ratios of the parameters jii are related to the structure of the Majorana neutrino mass 
matrix, 

= Q,(/,0)(/j0) . (2.21) 

An interesting example, which can account for the large mixing angles in the neutrino 
sector, is 'anarchy' [34] where Qj = 0(1). In this case one also has 

^ = 0(1) . (2.22) 

For the purpose of illustration we shall use in the following sections the 'democratic' case 
/^i = A*2 = /^3 where the parameters ^^i, ^zi and ^wi have no flavour dependence. Alterna- 
tively, one may consider 'semi- anarchy' (cf. [7]). Other examples can be found in [35]. 

In the following we shall consider gravitino masses between 100 GeV and 600 GeV, for 
which the assumed hierarchy m"^^ -C "^3/2 ^ "^Lft "^^^ only be a rough approximation. 
A more detailed treatment would have to incorporate mixings with heavy particles of the 
supersymmetric standard model. However, we find the operator analysis useful to illustrate 
the main qualitative features of gravitino decays. 
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3 Antimatter from gravitino decays 



The scenario of decaying gravitino dark matter provides, for a wide range of gravitino 
masses and lifetimes, a consistent thermal history of the Universe, incorporating successful 
primordial nucleosynthesis and successful baryogenesis through leptogenesis. Furthermore, 
if dark matter gravitinos decay at a sufficiently large rate, the decay products could be 
detected through an anomalous contribution to the high-energy cosmic ray fluxes. In this 
section we shall discuss the constraints on the gravitino parameters which follow from the 
observations of the positron fraction by HEAT and PAMELA and of the antiproton flux 
by BESS, IMAX and WiZard/CAPRICE. 

The rate of antimatter production per unit energy and unit volume at the position r 
with respect to the center of the Milky Way is given by 

Q(£,r-) = ^i(^§. (3.1) 

"^3/2T3/2 dE 

where dN/dE is the energy spectrum of antiparticles produced in the decay, which we 
calculated employing event generator PYTHIA 6.4 [36]. On the other hand, p(r) is the 
density profile of gravitinos in the Milky Way halo. For definiteness we shall adopt the 
spherically symmetric Navarro- Prenk- White halo density profile [37]: 

P^^) = / / / M2 ' (^■2) 

(r/rc)[l + (r/rj]^ 

with Po ~ 0.26 GeV/cm^ and — 20 kpc, although our conclusions are not very sensitive 
to the choice of the density profile. 

After being produced in the Milky Way halo, charged cosmic rays propagate in the 
Galaxy and its vicinity in a rather complicated way before reaching the Earth. Antimatter 
propagation in the Milky Way is commonly described by a stationary two-zone diffusion 
model with cylindrical boundary conditions [38]. Under this approximation, the number 
density of antiparticles per unit kinetic energy, f{T,f,t), satisfies the following transport 
equation, which applies for both positrons and antiprotons: 

= ^ = V • [K(T, r) V/] + ^ [b{T, r}f] - V • [^(0/] - 2M(^)r,,„/ + Q(T, r), (3.3) 

where reacceleration effects and non-annihilating interactions of antimatter in the Galactic 
disk have been neglected, since the primary particles produced by the gravitino decay 
rarely cross the disk before reaching the Earth. 

The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The 
diffusion coefficient K{T, r) is assumed to be constant throughout the diffusion zone and 
is parametrised by: 

K{T) ^KoP 7^^ (3.4) 
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where 13 = v/c and TZ is the rigidity of the particle, which is defined as the momentum 
in GeV per unit charge, TZ = p{GeV)/Z. The normahzation Kq and the spectral index 
5 of the diffusion coefficient are related to the properties of the interstellar medium and 
can be determined from the flux measurements of other cosmic ray species, mainly from 
the Boron-to-Carbon (B/C) ratio [39]. The second term accounts for energy losses due to 
inverse Compton scattering on starlight or the cosmic microwave background, synchrotron 
radiation and ionization. The third term is the convection term, which accounts for the 
drift of charged particles away from the disk induced by the Milky Way's Galactic wind. It 
has axial direction and is also assumed to be constant inside the diffusion region: V^(r) = 
Vc sign(z) k. The fourth term accounts for antimatter annihilation with rate Fann; when it 
interacts with ordinary matter in the Galactic disk, which is assumed to be an inflnitely 
thin disk with half-width h = 100 pc. Lastly, Q{T,r) is the source term of positrons or 
antiprotons, deflned in Eq. (13.11) . The boundary conditions for the transport equation, 
Eq. (13. 3p . require the solution f(T,r,t) to vanish at the boundary of the diffusion zone, 
which is approximated by a cylinder with half-height L = 1 — 15 kpc and radius R = 
20 kpc. A detailed study shows that there is some degeneracy among the parameters of 
the propagation model. This can be used to decrease the size of L and thereby reduce the 
antiproton flux from gravitino decays relative to the flux from spallation [40]. 

The solution of the transport equation at the Solar System, r = Tq, z = 0, can be 
formally expressed by the convolution 



where Tmax = "^3/2 for the case of the positrons and Tmax = "^3/2 — "^p for the antiprotons. 
The Green's function G{T, T') encodes all the information about astrophysics (such as the 
details of the halo proflle and the propagation of antiparticles in the Galaxy), while the 
remaining part depends on gravitino properties. Analytical and numerical expressions for 
the Green's function for the case of positrons and antiprotons can be found in [16]. 

Finally, the flux of primary antiparticles at the Solar System from gravitino decay is 
given by: 



where v is the velocity of the antimatter particle. 

The calculation of the high-energy cosmic ray fluxes from gravitino decay is hindered by 
a large number of uncertainties stemming both from astrophysics, encoded in the Green's 
function, and from particle physics. However, as we shall show below, present observations 
of the positron fraction and the antiproton flux constrain the parameters of the model well 
enough to make deflnite predictions on the diffuse gamma-ray flux. These predictions will 
be tested by the Fermi LAT results in the near future, thus providing a crucial test of the 
scenario of decaying gravitino dark matter. 

/^From the particle physics point of view, the solution to the transport equation depends 
on the following unknown quantities: the energy spectrum of positrons or antiprotons 




(3.5) 



V 



(3.6) 
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produced in the decay, the gravitino mass and the gravitino hfetime. The energy spectrum 
of positrons and antiprotons depends crucially on the i?-parity breaking interactions of the 
gravitino. It was shown in Section 2 that when the gravitino mass is large, the branching 
ratios for the dominant decay channels are predicted to be X]t-BR-(^3/2 ~^ hvi) ^ 1/4, 
X]jBR(V'3/2 — > Zvi) ~ 1/4, BR(?/'3/2 W^ef) ~ 1/2, while the branching ratio for 
"03/2 1^ is predicted to be much smaller. Therefore, in this limit the injection spectrum 
of antiprotons from gravitino decay is fairly model-independent, being just a function of the 
gravitino mass. This is not the case, however, for the energy spectrum of positrons, since 
the flavour composition of the final state depends on the flavour structure of the i?-parity 
breaking couplings, which cannot be predicted without invoking a model of flavour. 

The observation by PAMELA of an excess in the positron fraction at energies extending 
at least until 100 GeV implies a lower bound on the gravitino mass of ~ 200 GeV if this 
positron excess is interpreted in terms of gravitino decays. Besides, as shown in [10], there 
exists a theoretical upper bound on the gravitino mass of ~ 600 GeV in supergravity 
models with universal gaugino masses at the GUT scale, stemming from the requirement 
of successful thermal leptogenesis, yielding a relatively narrow range for the gravitino mass, 
777.3/2 — 200—600 GeV. If the positron excess observed by PAMELA is unrelated to gravitino 
decays, the gravitino mass may be as low as ~ 5 GeV without yielding overclosure of the 
Universe [7]. 

3.1 Constraints from the antiproton flux 

The interstellar antiproton flux from gravitino decay can be calculated from Eq. (13.51) using 
the corresponding Green's function. However, this is not the antiproton flux measured by 
antiproton experiments, which is affected at low energies by solar modulation. Under the 
force field approximation [41], the antiproton flux at the top of the Earth's atmosphere is 
related to the interstellar antiproton flux [42] by the simple relation 

2m,Tis + T,| J^^^^is)' ^3.7) 

where Tjs = Ttoa + 0f, with Tjs and Ttoa being the antiproton kinetic energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively, and ipp being 
the solar modulation parameter, which varies between 500 MV and 1.3 GV over the eleven- 
year solar cycle. Since experiments are usually undertaken near solar minimum activity, 
we shall choose ipp = 500 MV for our numerical analysis in order to compare our predicted 
flux with the collected data. 

As discussed in [16], the calculation of the antiproton flux from gravitino decay suffers 
from uncertaintes in the determination of the physical parameters in the propagation of 
charged cosmic rays in the diffusive halo, leading to uncertainties in the magnitude of fluxes 
as large as two orders of magnitude at the energies relevant for present antiproton experi- 
ments. The requirement that the total antiproton flux from gravitino decay be consistent 
with measurements gives a lower bound on the gravitino mass which strongly depends on 
the choice of the halo model. In the following we shall adopt the MED propagation model, 
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Figure 1: Antiproton fluxes for 7713/2 = 200, 400, 600 GeV in the MED set of propagation 
parameters that saturate the antiproton overproduction bound (see text). Dotted hues: 
antiproton flux from gravitino decays, dashed hues: secondary antiproton flux from spal- 
lation in the case of minimal nuclear cross sections, solid lines: total antiproton flux. The 
gravitino lifetimes are r3/2 = 7 x 10^^ s, 3 x 10^^ s and 1.5 x 10^^ s, respectively. 
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which provides the best fit to the B/C ratio and measurements of flux ratios of radioactive 
cosmic-ray species [39]. 

In order to determine the maximally allowed exotic contribution to the total antiproton 
flux, a precise knowledge of the secondary flux of antiprotons from spallation of high-energy 
cosmic rays on the protons and helium nuclei in the interstellar medium is necessary. 
Unfortunately, the determination of this secondary flux is also subject to uncertainties. 
First, the choice of the propagation model can change the prediction of the secondary flux 
by 10 — 20%. More importantly, the uncertainty in the nuclear cross sections for p-p, p-He, 
He-p and He-He collisions can change the prediction of the secondary flux by 22 — 25% 
above or below the central value [43]. 

A conservative upper bound on the antiproton flux from gravitinos is obtained by 
demanding that the total flux is not larger than the theoretical uncertainty band of the 
MED propagation model. This means that a 'minimal' dark matter lifetime for the MED 
model can be deflned by a scenario where the secondary antiproton flux from spallation 
is 25% smaller than the central value, due to a putative overestimation of the nuclear 
cross sections, and the total antiproton flux saturates the upper limit of the uncertainty 
band which stems from astrophysical uncertainties discussed above. This amounts to the 
requirement that the antiproton flux from gravitino decays should not exceed ~ 50% of 
the central value of secondary flux from spallation. 

Using the above prescription we flnd the following lower bounds for the gravitino life- 
time: 

T372"(200) ~ 7 X lO'S, r3";^"(400) ~ 3 X lO'S, r3%"(600) ~ 1.5 x lO'S, (3.8) 

where the numbers in parentheses correspond to the gravitino masses ^3/2 = 
200, 400 and 600 GeV, respectively. The corresponding antiproton fluxes from gravitino 
decay, the secondary antiproton flux from spallation and the total antiproton fluxes are 
shown in Figure [1] together with the experimental measurements by BESS, IMAX and 
WiZard/CAPRICE, and the uncertainty band from the nuclear cross sections in the MED 
propagation model. The minimal lifetimes (13.81) can be compared with the gravitino life- 
times needed to explain the PAMELA positron fraction excess, which will be discussed in 
the next section. 

3.2 Comparison with electron/positron fluxes 

Using the procedure described in the previous section it is straightforward to calculate 
the positron flux at Earth from gravitino decay, we shall adopt for deflniteness the MED 
propagation model [39], characterised hj 5 = 0.70, = 0.0112 kpc^/Myr, L = 4kpc, 
Vc = 12 km/s. Note that the sensitivity of the positron fraction to the propagation model is 
fairly mild at the energies where the excess is observed, since these positrons were produced 
within a few kiloparsecs from the Earth and barely suffered the effects of diffusion. 

To compare the predictions to the PAMELA results, we shall calculate the positron 
fraction, deflned as the ratio of the positron flux to the combined positron and electron 
flux, $e+/(^e+ + ^e-)- the backgrouud fluxes of primary and secondary electrons, as 
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Figure 2: Contribution from dark matter decay to the positron fraction and the total 
electron + positron flux, compared with data from PAMELA and HEAT, and ATIC, 
Fermi LAT and HESS, respectively; mg/a = 200 GeV, rg/a = 3.2 x lO^^ s, for W^l"^ decays 
pure electron flavour is assumed. The "Model 0" background is used, and for comparison 
with Fermi-LAT data 25% energy resolution is taken into account. 
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well as secondary positrons, we extract the fluxes from "Model 0" presented by the Fermi 
collaboration in [44], which fits well the low energy data points of the total electron plus 
positron flux and the positron fraction, and is similar to the MED model for energies above 
a few GeV [45] . Then, the positron fraction reads 

$|i^(T) + $;;+^(T) + $DM(T) + k $;;'!^(T) ' ^ ' ' 

where k = 0{1) is the normalization of the astrophysical contribution to the primary 
electron flux, which is chosen to provide a qualitatively good fit to the data. 

We now discuss the hypothesis that the PAMELA positron excess is due to gravitino 
dark matter decay. For this to be the case, the gravitino mass must be at least 200 GeV. 
As discussed above, the branching ratios into Standard Model particles will be essentially 
fixed for gravitino masses of a few hundred GeV. The decay W^i^ then has a 

branching ratio of ~ 50%, and the hard leptons that are directly produced in these decays 
may account for the rise in the positron fraction if a significant fraction of these leptons 
has electron or muon fiavour. 

Consider first the extreme case that the decays occur purely into electron fiavour. For 
^3/2 = 200 GeV, the PAMELA excess can then be explained for the gravitino lifetime 
r3^2(200) ~ 3.2 X 10^^ s, as illustrated by Figure [21 Note that this lifetime is a factor 2 
smaller than the minimum lifetime given in (13. 8p . which we obtained from the antiproton 
constraint. In other words, an interpretation of the PAMELA excess in terms of gravitino 
decays is incompatible with the MED set of propagation parameters once antiprotons are 
taken into account. Nevertheless, the MIN model and other sets of parameters that yield 
intermediate values for the antiproton fiux can easily be compatible with both the positron 
fraction and the antiproton-to-proton ratio observed by PAMELA. The situation is very 
similar for 7713/2 = 400 and 600 GeV. 

Figure [2] also shows the predicted total electron + positron fiux together with the results 
from Fermi and ATIC. Obviously, the "Model 0" presented by the Fermi collaboration 
in [44], cannot account for the present data, and the contribution from gravitino decays 
makes the discrepancy even worse. In particular, the data show no spectral feature expected 
for decaying dark matter. On the other hand, gravitino decays may very well be consistent 
with the measured total electron + positron fiux once the background is appropriately 
adjusted. This is evident from Figure [3] where the contribution from gravitino decays is 
shown in the theoretically well motivated case of fiavour democratic decays. The figure 
also illustrates that, depending on the gravitino mass, the dark matter contribution to the 
PAMELA excess can still be significant. 

An obvious possibility is that both, the total electron + positron fiux and the positron 
fraction, are dominated by astrophysical sources. For instance, for the gravitino mass 
^3/2 = 100 GeV we obtain from the antiproton fiux constraint r^™(100) ~ 1 x 10^"^ s. 
As Figure H] demonstrates, the contribution from gravitino decays to the total electron + 
positron fiux and positron fraction is indeed negligible. Nevertheless, as we shall see in the 
following section, the dark matter contribution to the gamma-ray fiux can be sizable. 
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Figure 3: Contribution from dark matter decay to the positron fraction and the total 
electron + positron flux, compared with data from PAMELA and HEAT, and ATIC, Fermi 
LAT and HESS, respectively; ^3/2 = 200, 400, 600 GeV with the minimal lifetimes fl3.8p 
(top to bottom); for W^l^ decays democratic flavour dependence is assumed. The "Model 
0" background is used, and for comparison with Fermi LAT data 25% energy resolution is 
taken into account. 
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Figure 4: Contribution of gravitino decays to positron fraction and total electron + positron 
flux for m3/2 = 100 GeV and r3/2 = 1 x 10 



27 



4 Predictions for the diffuse gamma-ray flux 

The total gamma-ray flux from gravitino dark matter decay receives two main contribu- 
tions. The first one stems from the decay of gravitinos in the Milky Way halo, 



dE 



2E^ dN^ 



halo "^3/2 dE 87rr3/2 Jios 



ph^\o{l)dl, (4.1) 



where dN^/dE is the gamma-ray spectrum produced in the gravitino decay. The integra- 
tion extends over the line of sight, so the halo contribution has an angular dependence on 
the direction of observation. 

In addition to the cosmological contribution, the total gamma-ray flux also receives 
a contribution from the decay of gravitinos at cosmological distances, giving rise to a 
perfectly isotropic extragalactic diffuse gamma-ray background. The flux received at the 
Earth with extragalactic origin is given by 



dE 



I dV MT. '^ ,r. _3 ' (4-2) 



eg "^3/2 Ji d{Ey) + QJQMy 
where y = I + z, z being the redshift, and 

Here, ^3/2, and f^A are the gravitino, matter and cosmological constant density pa- 
rameters, respectively, pc is the critical density and the present value of the Hubble 
parameter. 

The cosmological contribution is numerically smaller than the halo contribution. More- 
over, the flux of gamma-rays of cosmological origin is attenuated by the electron-positron 



14 




0.1 1 10 100 

E [GeV] 




0.1 1 10 100 

E [GeV] 




0.1 1 10 100 

E [GeV] 



Figure 5: Predicted gamma-ray flux for m^i2 = 200, 400, 600 GeV for tlie minimal 
lifetimes (13.81) . We assume decays purely into electron flavour here. We show both the 
background obtained by Sreekumar et al. as well as the background obtained by Strong, 
Moskalenko and Reimer. 
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pair production on the extragalactic background light emitted by galaxies in the ultravio- 
let, optical and infrared frequencies [46]. However, the flux of gamma-rays originating from 
the decay of dark matter particles in the halo is barely attenuated by pair production on 
the Galactic interstellar radiation field at energies below 10 TeV [47] . Thus, the total flux 
is dominated by the halo component, yielding a slightly anisotropic gamma-ray flux [13] 
which is compatible with the EGRET observations [14] . A detailed study of the prospects 
of detecting this anisotropy with the Fermi LAT is beyond the scope of this paper and will 
be presented in more generality elsewhere [48]. 

The gravitino decay produces a continuous spectrum of gamma-rays which is deter- 
mined by the fragmentation of the Higgs boson and the weak gauge bosons. On the other 
hand, there exists a gamma-ray line at the endpoint of the spectrum with an intensity 
which is model-dependent For our numerical analysis we shall use the typical branching 
ratio in this channel derived in Section 2, 



V 1^3/2 J 



for gravitino masses in the range from 100 — 600 GeV. In Figure [5] the predicted diffuse 
gamma-ray flux is shown for m3/2 = 200, 400, 600 GeV and the respective lower bounds 
(13.81) on the gravitino lifetime. These spectra correspond to upper bounds on the signal in 
gamma-rays that can be expected from gravitino dark matter decay. 

For comparison, we show two sets of data points obtained from the EGRET measure- 
ments of the diffuse extragalactic gamma-ray background using different models of the 
Galactic foreground emission. An analysis by Strong, Moskalenko and Reimer using a 
model, optimised to better simulate the Galactic diffuse emission, revealed a power law 
behaviour below 1 GeV, with an intriguing deviation from a power law above 1 GeV [17]. 
For our present analysis, we shall show both sets of results as the status of the extragalactic 
background is currently unclear. For the background obtained by Moskalenko, Strong and 
Reimer, the extragalactic component is described by the power law [17] 



7T-| / P \ -0.32 

6.8 X 10-^ ( ) (cm^ str s)-^ GeV. (4.4) 
\GeV / 



dE 



bg 

The earlier analysis by Sreekumar et al led to a less steep background [51], 

-0.1 



dE 



1.37 X 10-^ - — (cm^ str s)"^ GeV. (4.5) 
\GeV / 



^We neglect in our analysis the contribution to the gamma-ray flux from inverse Compton scattering of 
high energy electrons and positrons on the interstellar radiation field. We estimate that this contribution 
is peaked at energies smaller than 0.1 GeV and has an intensity ^ C'(10~^)(cm^ str s)~^ GeV [49], 

thus giving a negligible contribution to the total flux, which is constrained by EGRET to he ~ 
10~^(cm2 str s)"i GeV at = O.lGeV. This contribution, however, can be sizable for larger dark matter 
masses [49,50]. 
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Figure 6: Predicted gamma-ray flux for 7713/2 = 100 GeV and T3/2 = 1 x 10^ s. We show 
botli tlie background obtained by Sreekumar et al. as well as the background obtained by 
Strong, Moskalenko and Reimer. 

For comparison with the data points, we have averaged the slightly anistropic halo sig- 
nal over the whole sky, excluding a band of ±10° around the Galactic diskjl We have 
conservatively used an energy resolution a{E)/E = 15% as quoted by Fermi [52]. 

It is remarkable that for both choices of the extragalactic background, the antiproton 
constraint allows for a sizable deviation from a power law background if the gravitino mass 
is above 200 GeV. Therefore, if such a deviation with the proper angular dependence is 
observed by Fermi LAT, the scenario of gravitino dark matter will gain support. Further- 
more, the existence of a gamma-ray line at the end of the spectrum is predicted, with 
an intensity that, as discussed in Section 2, depends on the model of i?-parity breaking. 
This line could be observed by Fermi LAT in the diffuse gamma-ray background, but also 
by the ground-based Cherenkov telescopes MAGIC, HESS or VERITAS in galaxies such 
as M31 [13]. For smaller gravitino masses the gamma-ray line becomes more prominent 
whereas the contribution to the continuous part of the spectrum decreases. This is illus- 
trated in Figure 6 for m3/2 = 100 GeV. 

The observation of the discussed features in the diffuse gamma-ray spectrum might, if 
interpreted as the result of gravitino decay, open the exciting possibility of constraining 
the reheating temperature of the Universe. More concretely, the thermal relic abundance 

^ The halo signal would have a larger degree of anisotropy if the dark matter halo is not completely 
uniform but presents substructures, as suggested by N-body simulations of Milky- Way-size galaxies. In 
the present calculation we are interested in the average flux in the whole sky excluding the Galactic disk, 
which depends on the total amount of dark matter in this region and not on the way it is distributed. 
Therefore, for our purposes it is a good approximation to neglect substructures and to assume a smooth 
dark matter halo profile. 



17 



of gravitinos is given by [10,53,54] 



^'3/2' 



" -^■^[lO^^GeYj \ m3/2 y'viTeV/ ' ^ 



Therefore, imposing that the thermal abundance of gravitinos should not be larger that 
the total dark matter abundance, the measurement of the gravitino mass by Fermi LAT 
and the measurement of the gluino mass at the LHC imply the following upper bound on 
the reheating temperature of the Universe: 

which is saturated when all the dark matter gravitinos are of thermal origin. This 
bound has impoprtant implications for the scenario of thermal leptogenesis, which requires 
Tr ^ 10^ GeV, as well as for many inflationary scenarios. 



5 Conclusions 

In supersymmctric theories with small i?-parity breaking thermally produced gravitinos 
can account for the observed dark matter, consistent with leptogenesis and nucleosynthesis. 
Gravitino decays then contribute to antimatter cosmic rays as well as gamma-rays. We 
consider gravitino masses below 600 GeV, which are consistent with universal boundary 
conditions at the GUT scale. 

Gravitino decays into Standard Model particles can be studied in a model-independent 
way by means of an operator analysis. For sufficiently large gravitino masses the dimension- 
5 operator dominates. This means that the branching ratios into hu, Zv and W^U' are 
fixed, except for the dependence on lepton flavour. As a consequence, the gamma-ray flux 
is essentially determined once the antiproton flux is known. On the contrary, the positron 
flux is model-dependent. 

The gamma-ray line is controlled by the dimension-6 operator. Hence, it is suppressed 
compared to the continuous gamma-ray spectrum. Its strength is model-dependent and 
decreases with increasing gravitino mass. 

Electron and positron fluxes from gravitino decays, together with the standard GAL- 
PROP background, cannot account for both, the PAMELA positron fraction and the 
electron -|- positron flux measured by Fermi LAT. For gravitino dark matter, the observed 
fluxes require astrophysical sources. However, depending on the gravitino mass and the 
background, the dark matter contribution to the electron and positron fluxes can be non- 
negligable. 

Present data on charged cosmic rays allow for a sizable contribution of gravitino dark 
matter to the gamma-ray spectrum, in particular a line at an energy below 300 GeV. 
Non-observation of such a line would place an upper bound on the gravitino lifetime, and 
therefore on the strength of it!-parity breaking, restricting possible signatures at the LHG. 
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